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Abstract 

In order to model nonlinear viscous dissipative motions in solids, 
acoustical physicists usually add terms linear in E, the material time 
derivative of the Lagrangian strain tensor E, to the elastic stress ten- 
sor a derived from the expansion to the third- (sometimes fourth-) 
order of the strain energy density £ = £(tr E, tr E 2 , tr E 3 ). Here, it 
is shown that this practice, which has been widely used in the past 
three decades or so, is physically wrong for at least two reasons, and 
that it should be corrected. One reason is that the elastic stress tensor 
a is not symmetric while E is symmetric, so that motions for which 
er + cr T ^ will give rise to elastic stresses which have no viscous pen- 
dant. Another reason is that E is frame-invariant, while a is not, so 
that an observer transformation would alter the elastic part of the total 
stress differently than it would alter the dissipative part, thereby vio- 
lating the fundamental principle of material frame-indifference. These 
problems can have serious consequences for nonlinear shear wave prop- 
agation in soft solids, as seen here with an example of a kink in almost 
incompressible soft solids. 



1 Introduction and main statement 



Nonlinear elastic wave propagation is a subject of considerable interest for 
many scientific and industrial applications such as geophysical exploration, 
soft tissue acoustics, and the dynamics of rubbers, silicones, and gels. From 
a theoretical point of view the mathematics and mechanics of nonlinear wave 
phenomena is a classical yet still active subject of research where many out- 
standing problems are awaiting a definitive systematic treatment. 

In Physical Acoustics, the expansion of the strain energy density £ to 
include nonlinear corrections is often attributed to Landau. Hence at the 
"third-order" in the Lagrangian strain tensor E, we write 



where I k = tr E k (k = 1,2,3) are the strain invariants, A and \x are the 
(second-order) Lame coefficients, and A, B, C are the (third-order) Lan- 
dau coefficients. From an historical point of view, it is interesting to note 
that this notation by Landau seems to have appeared first in the Theory 
of Elasticity (1986), co-written with Lifshitz, as an unnumbered equation 
in an exercisepQ. In 1937, Landauj2] had already provided the third-order 
expansion of S, that time denoting the third-order elastic constants by the 
letters P', Q', R'. That same year, Murnaghan[3j also proposed a third-order 
expansion using different strain invariants, which is still in use today (mostly 
by geophysicists). However the paternity of the third-order expansion can be 
traced further back in time, at least to a 1925 paper by Brillouin[4J, who in 
fact seems to have also been the first to use the letters A, B, C for third-order 
elastic constants. 

In any event, Landau and Lifshitz denote by Xk the Lagrangian coordinate 
components, by Uk the mechanical displacement components, and they derive 
the equations of motion in Cartesian coordinates as follows, 



\AI Z + Bhh + \cr 
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(1) 



do-jk 
dx k 




(2) 



Po 



dt 2 ' 



where po is the mass density in the undeformed configuration, and the stress 
tensor is related to the strain through 



(It is easy to recognize that this is the "first Piola-Kirchhoff stress tensor" 
of Classical Continuum Mechanics.) For the reader's convenience, we recall 
that the full Lagrangian strain tensor E has Cartesian components 

1 / dui du k duj duj 



2 \dx k dxi dxidx k J ' 

and can be rewritten as 

Eik = \ (FjiFjk — Sik) ■ (5) 

Here, F ik are the Cartesian components of F, the deformation gradient, 
defined as 

F ik = t— = o ik + — , (6) 
ox k dx k 

where yi(x k ) are the space coordinates of the current position of material 
point x k , for the deformation given by yi(x k ) = Xi + Ui(x k ). 

In another part of their book (more precisely: at the end of §34 in Chapter 
5) Landau and Lifchitz model viscous dissipation by adding to the elastic 
stress tensor cr a "viscosity stress tensor" cr', with Cartesian components 



a' ik = 2 V (E ik - U ik E u ) + (E u 5i 



(C-iv)Eu5 ik + 2rjE, 



(7) 



ik 



where rj > and £ > are the shear and bulk viscosity coefficients, re- 
spectively, and the superposed dot denotes the time derivative. Note that it 
is completely unambiguous from the context of this part of their book that 
Landau and Lifshitz are speaking here of the infinitesimal theory of visco- 
elasticity, and that E in Eq. ([7]) is the time derivative of the infinitesimal 
strain tensor (i.e. the linear part of Eq. Q), see their unnumbered equation 
between (34.2) and (34.3). It might seem at first glance that adding cr' to cr, 
and taking the strain components to be those of the finite Lagrangian strain 
Q instead of the infinitesimal strain, would be a first, logical step towards 
the inclusion of nonlinear dissipative effects. However there are at least two 
problems with this seemingly anodyne approach. 

One problem is that cr' is symmetric while cr is not. Then the resultant 
total stress tensor cr + cr' has a visco-elastic symmetric part: (cr4-cr T )/24-cr', 
but a purely elastic anti-symmetric part: (cr — cr T )/2, where the superscript T 
denotes the transpose. It is thus impossible to model antisymmetric viscous 
stress effects in motions and boundary conditions with this formulation. 
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The other problem is that cr and a' are made objective, or frame- 
indifferent, in two different and irreconcilable ways. Indeed, it is well known 
that there are transformation rules to follow in order to ensure that the direc- 
tions associated with a tensor are unaltered by an observer transformation. 
To summarize, two observers, one associated with a frame with position x 
and time t and the other associated with a frame with position x* and time t* 
are said to be equivalent when they are connected by (see e.g. Chadwick[5]), 

x* = c(t) + Q(t)x, t* =t-a, (8) 

where Q is a proper orthogonal tensor, c is a vector, and a is a constant 
scalar. The transformation rule for the elastic part of the stress tensor ^ is 
that cr transforms into cr given by 

a = Q(t. (9) 

in order to ensure frame-indifference, whilst for the viscous part (JTJ) , it is 

c7 = cr', (10) 

because E - and hence E — is observer-invariant (this is not so obvious 
from Eq. (|4]) , but becomes more so in view of Eq. ^ , as shown in Ref . [5] ) . 
Clearly, it is not possible for the composite stress tensor cr + cr' to comply 
with both requirements simultaneously and appear the same to two equiv- 
alent observers. This state of affair violates the fundamental principle of 
objectivity, a cornerstone of rational mechanics, e.g. see Gurtin[6]. 

Now, there are two ways to reconcile cr and cr' in order to construct a 
stress tensor coherent with respect to symmetry and objectivity. One way 
is to make cr behave like cr' , by replacing it with F~ l cr, the (symmetric) 
second Piola-Kirchhoff stress tensor; the other way is to make cr' behave like 
cr, by replacing it with cr' = Fcr' . 

Evidently, both courses of action are eventually equivalent. In weakly 
nonlinear elasticity theory, governing equations have been around for cr longer 
than for cr' and we henceforth take advantage of these and concentrate on the 
consequences of modifying cr' rather than cr. We thus propose the following 
form for the viscous stress, cr' = Fcr', with components 

®ik = Fij a jk 

= (C - | rj) F lk E n + 2nF tj E jk ^ 

= (c-fi)(*.+g^(«. + g)^. 
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and investigate the consequences of using cr' instead of a' in the equations 
of motion. 

In passing, it is interesting to show the form taken by the constitutive 
equation for the Cauchy stress tensor r', related to cr' by 

tfk = v'isFks, (12) 
where j = det F. For this we recall that [5] 

Eik = FudijFjk, (13) 
where dij is the Cartesian component of d, the Eulerian stretching tensor, 

1 / dvi dvj\ ,„ t . 

d » = 2 {4 + at) ■ (14) 

and that the (symmetric) left Cauchy-Green strain tensor b is defined as 

bji = FjiFu, (15) 



so that, in view of Eq. (13), the expression for En can be written as 



En = FudijFji = bjidij = tr(bd). (16) 

Thus, Eq. (11) yields 

°ik = (C-h)FikE u + 2riF ij E jk 

= (C - lv) F ik bfldn + 2-qFy F tj d ls F sk (17) 
= (C - Irj) F ik + 2r)bud la F ak . 



From Eqs. (12), (13) and (|16j) we compute the viscous part f[ k of the Cauchy 
stress as ^ f 

= (C - fry) F ik F mk bijdji + 2rjbud ls F sk F mk (18) 

= (C - lv) hmhidis + 2r)hidi s b sm , 
or, in absolute notation, 

jt' = (C - Irj) btr(bd) + 2r]bdb. (19) 

It is quite clear that this is a polynomial isotropic function of the two tensor 
variables b and d, which can be shown to be a special case of the universal 
representation for such functions (see, e.g., the classical treatise by Truesdell 
and Noll [TJ Sect. 13] and the paper by Rivlin|8j). It is reassuring to notice 



that expression (19) for the viscous part of Cauchy stress f' falls perfectly 
within the class of constitutive equations for finite viscoelasticity. This makes 
much more plausible our proposal (11) for a modified version of Eq. 0. 



2 Wave motion 



When it comes to study wave propagation, attention is usually focused on 
special motions, for which the general equations of motion simplify greatly. 
One might then be led to believe that no relevant corrections should be made 
on the governing equations as a consequence of the substitution of Eq. Q 
with the more appropriate version Eq. (11). In this section we show that 



this is not the case, indeed, and that some relevant and quite complicated 
changes should be made to the propagation equations, even for simple (bulk) 
motions. 

We discuss the propagation in a soft material of a transverse wave de- 
scribed by a displacement field Ui(xk,t) of the form 



u\ = u(z, t), 



u 2 = u 3 = 0, 



(20) 



where z is the third Cartesian coordinate. This is the choice made by Zabolot- 
skaya et al. [9] , who added a viscous part in the form ([7| to an elastic consti- 
tutive equation (see their Eq. 34). 

Straightforward computations show that 




1 



E 















<r 



(21) 



(here, for compactness, we denote partial derivatives with subscripts, while 
in other formulas we use the more explicit notation, so that a comparison 
with some relevant references is made easier). 

The time derivative of E and its trace are thus given by 



E 



and, moreover, 



FE 






u zt 









U z t 



2u z u zt 



tr(E) 



U z Uzt, 



(22) 



U z 


1 






Uzt 



U zt 



2u z u z t 



u z u z t u zt + 2u 2 z u zt 




Uzt 



zt 



(23) 
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Now, since by Eq. (11) 



(C- h)ti(E)F + 2r]FE, 



(24) 



we find that the Cartesian components of the viscous part cr' of the Piola- 
Kirchhoff stress tensor are given by 



(C + \rj)u z u zt 

(C 

Wzt 



r]u zt + (C + \r})u\u zt 

\f])u z u zt 

(C + \rj)u z u zt 







(25) 



In their Section III, Zabolotskaya et al.[9j investigate the propagation of 



waves described by Eq. (20) in an almost incompressible soft material with 
elastic strain energy 



S = /2l 2 + \Ah + Dll 



(26) 



where //, A, D are second-, third-, and fourth-order elastic constants, respec- 
tively. Terms of order higher than the third in the strain du/dz are then 
neglected and the governing equations reduce to the single equation 



Po 



d 2 u 



dt 2 



d 2 u ^ d fdu\ 
^ dz 2 ^ dz \dz J 



(27) 



where po is the constant density and 7 = fi + \A + D, see Ref. [9J Eq. 12]. 

Next, in Section V of Ref. [9], a viscous stress is added, in the form of a' ik 
as defined here by Eq. ([7]) , and an appropriate modification for the governing 
equation (27) is deduced, which leads to the introduction of an additional 



viscous term: 



d 2 u 
di 2 



d 2 u 







P 



dz 2 ^ dz 



du 

dz 



+ 7] 



0* 
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dz 2 dt 



(28) 



elastic 



This equation can also be found in the Refs. jTOJ, [TTJ, [12] , for instance. Our aim 
here is to verify which modifications would be required for the wave equation 



(27) if the added viscous Piola-Kirchhoff stress were defined by Eq. (11), as 
we suggest, rather than by Eq. ([7]). In other words: what difference does it 
make to the wave propagation equations to consider such a stress tensor cr' ik , 
rather than a' ik l 
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In order to answer this question we first compute the divergence of cr', as 



(Divff')i 
(Div<7-') 2 



d 3 u 4 d 



d ( du d 2 u 



( d 2 u \ ( du 
\dzdi) [dz 



(29) 



(Div^) 8 = (C+f^( 



dz dzdt) 



Here we remark that the first and last of these components lead to two 
non-trivial equations of motion, in contrast to the situation encountered in 
Ref. [9], where there was only one equation. The second equation here 
may however be made to disappear if we were to consider soft solids to be 
perfectly incompressible, and had thus to introduce an arbitrary Lagrange 
multiplier (for instance, see Ref. [13] for a rational inclusion of perfect in- 
compressibility in weakly nonlinear elasticity, and Ref. [H] for a derivation 
of the equations of motion in compressible and incompressible materials and 
the possible decoupling of longitudinal from transverse waves.) 
Next, we readily derive the wave equation 



d 2 u 



Po 



d 2 u 

dz A 



+ 7 



d_ 

dz 



du 

dz 



+ 7] 



<9 3 



dz 2 dt 



+ (c + h) 



d_ 

dz 



f d 2 u \ (duX 
[dzdt J \dz) 



elastic 



v 

viscous 



(30) 

The viscous part of this equation is remarkably more complex than the sim- 



pler term rjd 3 u/dz 2 dt shown in Eq. (28), which is Eq. (36) of Ref. [9]. Notice 
that the last addendum in Eq. (|30[) is of third order in the strain du/dz 



exactly as the nonlinear elastic term d/dz (du/dz) . Thus, since Eq. (27) 



was obtained in Ref.[9j by neglecting terms of order higher than three, it is 



our opinion that, in principle, the additional viscous term in (30) should be 



kept and considered for a coherent discussion of wave propagation. Dimen- 
sional analysis reveals that this term would be negligible if fifth-order elastic 
constants were much larger than lower-order constants. Although this might 
be the case for some solids, it can be shown by following the steps presented 
in Ref. [13] and Ref. [18] that for almost incompressible solids, all elastic con- 
stants are of the same order as /i (see Ref. [19] for experimental evidence). 



To illustrate the potential influence of the extra term in Eq.(30), we 
focus on a staple of acoustic nonlinearity: the finite amplitude, traveling 
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kink solution. First we rewrite the equation of motion for the strain 
w = du/dz, as 

d 2 w d 2 w . d 2 3 d 3 w , 4 v d 2 ( 2 dw\ 



Then we perform the following changes of variables and of function: 

{ = ^«, r = ^ W"' = ./^, (32) 

T] T] \ fi 

to obtain 

d 2 W d 2 W d 2 ,„ r , 3 d 3 W n (h , , d 2 ( TTr2 dW\ . . 

Looking for a traveling wave solution in the form 

W{£, r) = Vc 2 - 1 where x = (c 2 - l)(r - £/c), (34) 

and c > 1 is the arbitrary, non-dimensional speed, we find the following 
equation for the amplitude to, 

J' = ( w 3) " + w /» + Q « ^ where a = (/i/7) ^ + _ x y (35) 

It can be integrated twice for a finite kink with tails such that u(— oo) = 0, 
w(oo) = 1, u/(±oo) = u/'(±oo) = 0, to give 

(l -cj 2 ) = (i + a uj 2 )uj'. (36) 



Now the variables of this equation can be separated, and by centring the kink 
so that u(0) = 1/2, we obtain its inverse definition, 

* = -lni-[f(l-u, 2 )]^. (37) 

It is a simple matter to construct the to — x curves and then to invert them 
to generate the x — to curves. Explicit inversions include 



at a = 0, 1, respectively. At a = 0, there is no extra term in Eq.(30); as soon 
as a 0, the difference between the incorrect formulation and the proper 
formulation of viscous effects is felt, as illustrated by Figure 1. 
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Figure 1: Finite amplitude travelling transverse strain kink in a viscous soft 
solid. The equation of motion has been fully non-dimensionalized. The pa- 
rameter a is a measure of the consequence of properly incorporating viscous 
effects in the formulation of the equations of motion. At a = (dotted 
curve), the viscous formulation does not allow for the principle of objectiv- 
ity to be respected. As a grows (here, a = 1,2,3 in turn), the wave front 
displays a gentler slope and the wave is more attenuated. 



3 Concluding remarks 



Viscous stresses are often introduced in the formulation of nonlinear wave 
motion problems in order to prevent the formation of shocks, because they 
confer a parabolic character to the equations of motion for continua. How- 
ever, as pointed out earlier by Antman [16], the choices made historically in 
the literature for these stresses sometimes turn out to be physically unac- 
ceptable because their material responses are affected by rigid motions. Here 
we showed that simply extending the linear Kelvin- Voigt model of differen- 
tial visco-elasticity from linearized to finite elasticity is not a straightforward 
process. In particular, one choice is to take the strain-rate effects to be de- 
scribed by the time-derivative of the full Lagrangian strain instead of the 
infinitesimal strain. With that choice must come great care in formulating 
a corresponding viscous stress which obeys the fundamental principle of ob- 
jectivity in the same manner that the elastic stress does. We saw here that 
this compatibility can be achieved by pre-multiplying the linear expansion 
of the viscous stress in terms of ii by the deformation gradient F, which 
complicates greatly the equations of motion, as illustrated in Section [2j 

An alternative constitutive assumption for modeling viscous effects is to 
take the Cauchy stress tensor to be linear in d, the Eulerian stretching tensor 



defined in Eq. (14). That assumption is perfectly coherent with the funda- 



mental principles of mechanics, including frame-invariance, and is aligned 
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with modeling of viscous effects in Fluid Mechanics and the emergence of 
the Navier-Stokes equations. It has also been used to model nonlinear wave 
propagation in solids[Hl [El HZj- It does not complicate the equations of 
motion excessively. 
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